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ABSTRACT 

The Herschel ATLAS is the largest open-time key project that wiU be carried out on the Hcrschel 
Space Observatory. It wiU survey 510 square degrees of the extragalactic sky, four times larger than 
all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the 
survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, 
and the six major science programmes we are undertaking. Using new models based on a previous 
submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in 
other wavebands. 
Subject headings: 
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1. INTRODUCTION 

Approximately half the energy emitted since the 
big bang by all the objects in the Universe has 
been absorbed b y dust and then reradiated between 
60 and 500 urn ^Dwek et al. 1998t iFixsen et al. 19981 : 
[Driver et al. 20081 ) , a wavelength range in which the Uni- 
verse is still largely unexplored (Fig. 1). On the short- 
wavelength side of this waveband, the whole sky was sur- 
veyed at 60 and 100 /im by IRAS in the 1980s. However, 
almost all of the tens of thousands of galaxies detected by 
IRAS were spirals and starbursts in the nearby Universe 
(z < 0.1), and IRAS revealed little about the dust in 
other gal axy popula tions, especially early-type galaxies 
(jBregman et al. 199 8'). Even in the late-type galaxies, 
only the small fraction of the the dust warm enough to 
radiate significantly in the far-infrared was detected by 
IRAS. Devereux and Young (1990), for example, showed 
that the gas-to-dust ratio estimated from IRAS measure- 
ments alone was ~10 times greater than the standard 
Galactic value, implying that ~90% of the dust in galax- 
ies was effectively missed by IRAS. ISO and Spitzer, with 
their long wavelength (170 /im) band, suffered less from 
this problem b ut will still have effe ctively missed any dust 
with T < 15K (jBendo et al. 20031 ). Notwithstanding the 
successes of IRAS, ISO and Spitzer, most of the wave- 
band from 60 to 500 fim is still virtually terra incognita, 
and the only survey of a large area of the extragalactic 
sky at a wavelength beyond 200 /xm is the one recently 
carried out by the Herschel pathfinder experiment, the 
Balloon Large Area Survey Telescope (BLAST), which 
covers ~ 20 deg^ (jPevlin et al. 2 009). 

The lack of a survey covering a large area of sky 
in the submillimetre waveband (in this paper defined 
as 100 /im < A < 1 mm) has left us in some ways 
knowing more about dust in the distant, early Universe 
than in the Universe today. The surveys that have 
been carried out in the submillimetre waveband with 
ground-based telescopes — at 450 and 850 fim with 
the SCUBA camera on the J ame s Clerk Maxwell 
Telescope (iHughes et al. 19981: lEales et al. 19991: 
ICoppin et al. 2006f) . at 1.2 mm with MAMBO 
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Fig. 1. — The extraglactic b ackground radiation as a function of 
wavelength IIDole et al. 20061 ) 

on the IRAM 3 0- m telescope (iGreve et al. 20041: 
iBertoldi et al. 20071 : IGreve et al. 20081 ) and at 1.1 
mm wit h AZTEC on the Jan ies Clerk Maxwe ll Tele- 
scope (jPerera et al. 20081: iScott et al. 2008D — have 
been of very small areas of sky, covering ~ 1 deg^ of 
sky in total. Because of the unusual submillimetre 
'K-correction'^^, these surveys have mostly detected 
sources at very high redshifts (z >z 1). They have led 
to the important discovery that there is a population of 
luminous dust-ens hr ouded galaxies in th e early Universe 
(jSmail et al. 1997t iHughes et al. 1998f) . which many 
author s have suggested are the ancestors o f ellipticals 
today (jScott et al. 20021 : iDunne et al. 20031 ) . but they 
have told us relatively little about the evolution of 
the Universe since z ~ 1 (the last 8 billion years) and 
almost nothing at all about dusty galaxies in the nearby 
Universe. 

Two basic things one would like to know about the 
nearby Universe are the submillimetre luminosity func- 
tion and the dust-mass function(thc space-density of 
galaxies as a function of dust mass). These functions 
are important for many reasons, includ ing tests of semi- 
analytical models of galaxy formation (|Cole et al. 20001 : 
iBaugh et al. 200"5l ) and, by comparison with the same 
functions at high redshift, accurate measurements of the 
amount of evoluti on that is oc curring in the submil- 
limetre waveband (jPunne et al. 2003) . Unfortunately, 
our knowledge of these functions is still extremely poor 
because of the limitations in areal coverage and sen- 
sitivity of previous submillimetre telescopes. Until re- 
cently the only estimates of the local submillimetre lu- 
minosity function were one base d on 55 galaxies de- 
tected in an ISO 170-/im survey (Takeu chi et al. 20061 
and ones based on SCUBA 850-/ini o bservations of ^200 
galaxies selected in o ther wavebands ()Dunne et al. 2000t 
iVlahakis et al. 20051 ). There are now direct estimates 
from the BLAST results of the local luminos ity function 
at 250, 350 and 500 /xm (jEales et al. 200 9t). but their 
accuracy is limited by the small number of low-redshift 
sources detected in the BLAST survey: ~ 30 at z < 0.2. 

Beyond a redshift of ~ 1, as the redshift increases, the typi- 
cal spectral energy distribution of a dusty galaxy means that the 
effect of increasing luminosity-distance on the brightness of the 
galaxy is compensated for by the increasing rest-frame luminosity 
of the galaxy. The consequence is that the galaxy's flux density is 
approximately independent of redshift. 
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The lack of a large-area survey capable of measuring 
the dust content and dust-obscured star formation in 
large numbers of galaxies in the local Universe has been 
especially galling for submillimetre astronomers in the 
light of the success of their colleagues working in other 
wavebands. The Sloan Digital Sky Survey and the 2dF 
Galaxy Redshift Survey have led to a revolution in our 
understanding of the distribution of galaxies in the local 
Universe, and the relationships between their present 
star-formation rate, star-formation history, stellar 
ma ss, morphology and environment (Lew is ct al. 2 002; 

Kauffmann et al. 2003al: |K auffma nn et al. 2d03bl : 

Heavens et al. 20041: iBalogh et al. 2004. ). However, all 
the studies that have used these impressive datasets 
to investigate the physics and ecology of the galaxy 
population have been forced to ignore the dust phase 
of the interstellar medium and star formation that is 
heavily obscured by dust, because the IRAS survey was 
only sensitive enough to detect a small percentage of the 
galaxies in the redshift surveys (1.8 per cent of the SDSS 
galaxies - Obric et al. 2006), and it missed 90% of the 
dust in the detected galaxies because of its insensitivity 
to cold dust (Dcvcrcux and Young 1990). 

The launch of the Herschel Space Observatory, which 
occurred on May 14th 2009, has the potential to dramat- 
ically increase our knowledge of dust and dust-obscured 
star formation, especially in the nearby Universe. Her- 
schel has two main cameras: SPIRE, which will be able 
to image the sky simultaneously at 250, 350 and 500 
fim ( Griffin et al. 2007. ) . and PACS, which will be able 
to image the sky in two bands si multaneously, either 7 
and 170 ^j,m or 110 and 170 /xm (jPoglitsch et al. 20061 ). 
Herschel will have much better angular resolution (~18 
arcsec at 250 /im) and sensitivity than previous observa- 
tories, and the spectral coverage of SPIRE will make it 
possible to carry out the first large-area surveys in this 
virtually unexplored part of the electromagnetic spec- 
trimi. The SPIRE bands will also make it possible to 
detect the cold dust that was missed by earlier obser- 
vatories. Herschel is also better suited for investigating 
the local Universe than the submillimetre surveys that 
will soon be carried out from the ground, in particular 
the SCUBA-2 and LABOCA surveys, because these will 
operate mostly at 850 /um, where low-redshift galaxies 
are intrinsically faint, whereas the Herschel bands span 
the peak of the typical spectral energy distribution of a 
galaxy in the nearby Universe. 

In this paper, we describe the largest project that will 
be carried out with the Herschel Space Observatory in 
'Open Time', the time available for competition within 
the international astronomical community'^^. The Her- 
schel Astrohysical Terahertz Large Area Survey (the Her- 
schel ATLAS or H- ATLAS) will be a survey of 510 deg^ 
of sky in five photometric bands. This is ~8 times larger 
than the coverage of the next largest (in area) Herschel 
extragalactic survey, HERMES (Oliver et al. 2009). The 
main scientific goal of the H- ATLAS is to provide mea- 
surements of the dust masses and dust-obscured star for- 
mation for tens of thousands of nearby galaxies, the far- 
IR/submillimetre equivalent to the SDSS photometric 

This forms about 2/3 of the total observing time on Herschel, 
with the remaining 1/3 being time reserved for the teams that built 
the instruments — 'Guaranteed Time'. 



TABLE 1 
H-ATLAS Fields 



Name 


Centre 


RA width" 


Dec width" 


NGP 


13 18 00, 29 00 00 


15 


10 


GAMA A 


09 00 00, 00 00 00 


12 


3 


GAMA B 


12 00 00, 00 00 00 


12 


3 


GAMA C 


14 30 00, 00 00 00 


12 


3 


SGP A 


02 26 48, -33 00 00 


11 


6 


SGP B 


23 15 36, -32 54 00 


31 


6 


Note. — 


Reading from the left, the eolumns are: 


the name of the 



field; the central position of the field; the width of the field in degrees 
in RA; the width of the field in degrees in declination. 
^ The precise coverage of the fields and their orientation on the sky will 
depend on exactly when they are observed. See §6 for more details. 

survey. However, the H-ATLAS has many other science 
goals ranging from the investigation of the point sources 
that will be detected by the Planck Surveyor to a study 
of high-latitude galactic dust. 

The arrangement of this paper is as follows. In Section 
2 we describe the basic parameters of the survey. In Sec- 
tion 3 we present predictions of the number and redshift 
distribution of the sources that will be detected by the 
H-ATLAS. In Section 4, we describe the complementary 
data that exists or will soon exist for the H-ATLAS fields. 
In Section 5 we describe the six main H-ATLAS science 
programmes. In Section 6 we describe the detailed sur- 
vey strategy, including issues that will be of interest to 
the general community, such as our plans for the release 
of data products. We everywhere assume the cosmologi- 
cal parameters for a 'concordance universe': Qm — 0.267 
and flA = 0.762. 

2. THE BASIC PARAMETERS OF THE SURVEY 

The H-ATLAS has been allocated 600 hours of time, 
making it the largest key project that will be carried 
out with Herschel in Open Time. For all of our science 
goals (§5), the final sensitivity of the survey is not criti- 
cal and it is more important to survey the greatest pos- 
sible area of sky. We have therefore chosen to use the 
maximum possible scan rate for the telescope (60 arc- 
sec sec~^). For our first science programme (§5.1), it is 
important to make observations with PACS and SPIRE, 
and therefore we have chosen to use the Herschel ob- 
serving mode that allows simultaneous observations with 
the two cameras: Parallel Mode (PMode). Of the two 
possible combinations of photometric bands for PACS 
fPoglitsch et al. 2006.) . we have chosen to observe at 110 
and 170 /xm rather than at 70 and 170 fim mostly on 
the grounds of sensitivity; the noise at 70 and 110 ^m 
should be fairly similar but galaxies, even at low red- 
shift, are generally brighter at the longer wavelength. 
Although this combination will be worse for estimating 
the temperature of the dust, our models suggest that we 
will still be able to obtain useful measurements of the 
temperature of the dust in low-redshift galaxies. 

With an eye on the legacy value of the H-ATLAS, we 
have chosen to observe fields in the northern and south- 
ern hemispheres and on the celestial equator. Other than 
that, we have chosen our fields to maximise the amount 
of complementary data and to minimize the amount of 
confusing emission from dust in the Galaxy, this last de- 
termined from the IRAS 100 /Ltm maps. The fields, which 
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Fig. 2. — The positions of the ATLAS field, shown as white blocks, superimposed on tlie IRAS 100 fim map of the sky, which traces the 
distribution of galactic dust. Figure 2(a) shows the northern galactic cap and Figure 2(b) shows the southern galactic cap. The colour 
coding for the lines is as follow: solid green lines — RA and dec; dotted green lines — ecliptic latitude and longitude; cyan — KIDS/VIKING 
area; yellow — SDSS area; blue — 2dFGRS area; magenta — area of the Dark Energy Survey; magenta/blue dashed — area covered by the 
South Pole Telescope. 



are shown in Figure 2 and listed in Table 1, are: 

• One field close to the northern galactic pole with 
an area of 150 deg^ (henceforth the NGP field) 

• Three fields, each of 36 deg^ in area, coinciding 
with the fields being surveyed in the Galaxy And 
Mass Assembly redshift survey (Driver et al. 2009) 
(henceforth the GAMA fields) 

• Two fields with a total area of 250 deg^ close to the 
south galactic pole (henceforth the SGP fields) 

The total survey covers 510 deg^. The angular resolu- 
tion (full-width half maximum) of the observations will 
be approximately 8, 12, 18, 25 and 36 arcsec at 70, 110, 
250, 350 and 500 fim, respectively. The 5f7 sensitivi- 
ties that should be reached in the five bands are 67 mJy 
at 110 /im, 94 mJy at 170 /xm, 45 mJy at 250 /im, 62 
mJy at 350 /Ltm and 53 mJy at 500 /xm. These sensi- 
tivities have been estimated using the current version of 
the Herschel observation planning package HSpot. We 
have assumed that the sources will be unresolved by the 
telescope beam, which will be true for most sources but 
will not be true for the closest sources. We have also 
neglected the effects of emission from dust in the Galaxy 
and of confusion from extragalactic sources. The im- 
portance of these is still uncertain, but the importance 
of confusion s hould be l ess than for the deeper Herschel 
surveys ((Oliver et al. 20 09). 

3. PREDICTIONS 

The number of sources that will be found in any of the 
Herschel surveys is uncertain because they cover regimes 
of flux density and wavelength that have never been ex- 
plored before. There is a particular problem in the wave- 
length range 200 < A < 500 /zm because no survey, un- 
til the very re cent survey with t he balloon-borne tele- 
scope BLAST (iPe vlin et al. 2009 f). has been carried out 
at these wavelengths. Models based on surveys at shorter 



wavelengths, such as IRAS, are likely to underestimate 
the number of sources in the submillimetre band because 
these surveys miss the col d dust tha t radiates strongly 
at the longer wavelengths (jPunne fc Bales 2001). 

We have used two models to predict the proper- 
ties of the sources that will be found in the H- 
ATLAS. First, the model of Lagache and collaborators 
(jLagache et al. 20031 : iLagache et al. 20041 ) IS an empir- 
ical evolution model in which an analytical form is as- 
sumed for the evolution, and the parameters of the 
model are adjusted to fit the observational data in the 
far- infrared and submillimetre wavebands, in particu- 
lar the source counts at 24, 70, 160 and 850 /im and 
the spectrum of the cosmic background radiation. Sec- 
ond, we have developed a new empirical evolution model 
based on the SCUBA Local Universe and Galaxy Survey 
punne et al. 2000^ 'Vlahak is et al. 20 05). Both models 
give ade quate agreem ent with the results of the BLAST 
survey (|Dve et al. 2009 !). although the comparison be- 
tween the results and the models is still at a preliminary 
stage because of th e large amount of source confusion in 
the BLAST survey (|Eales et al. 2009( 1. We now describe 
the new model in more detail. 

The model is based on the sample of 104 galaxies ob- 
served by Dunne et al. (2000) with the SCUBA sub- 
millimetre camera as part of the SCUBA Local Uni- 
verse and Galaxy Survey (SLUGS). These galaxies form 
a statistically- complete sample ab ove a flux limit at 60 
/im of 5.24 Jy (! Dunne et al. 2000l ). All of these galaxies 
were detected at 850 /tm and many at 450 /tm. Dunne 
and Eales (2001) present a simple two-component dust 
model that fits the 60, 100, 450 and 850 nm flux mea- 
surements. For the galaxies without measurements at 
450 //m, it is possible to determine the parameters of 
the model by making the additional assumption that the 
temperature of the cold dust component is 20 K, which 
is the average of the estimates for the galaxies that do 
have complete flux measurements. The SLUGS sample is 
still the only large sample of galaxies for which there are 
empirical spectral energy distributions that extend from 
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Fig. 3a. — On the left is the redshift distributions predicted at 110 fj,m by the models of Lagache et al. (2003,2004); on the right is the 
redshift distribution predicted at the same wavelength by the model described in the text. 
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Fig. 3b. — Same as in (a) except at 170 /Lim. 
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Redshlft Redshift 

Fig. 3c. — Same as in (a) except at 250 //m. 



1 2 
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Fig. 3d. — Same as in (d) except at 350 jum. 
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Fig. 3e. — Same as in (e) except at 500 /im. 



the far-infrared to submiUimetre waveband. More sophis- 
ticated attempts to use the SLUGS data to predict the 
local luminosity function suggest that the simple method 
we present here leads to an underestimate of the lo- 
cal luminosity function (jVlahakis at al. 20051 ). However, 
the local submiUimetre luminosity function is sufficiently 
uncertain that we are not too concerned about this — 
determining the luminosity function is one of the goals 
of the H- ATLAS — and the simple modelling method we 
present below has the advantage that it is possible to 
make predictions for the optical, radio and other proper- 
ties of the galaxies we will detect with H- ATLAS. 

We can use the SLUGS sample to predict the source 
counts in any submiUimetre waveband in a straight- 
forward way. Let us, for example, predict the num- 
ber of sources in the SPIRE 250 fim band. If we as- 
sume 'number-density evolution', in which the number 
of sources of a given luminosity changes with redshift, 
the number of sources above a given 250 fim flux density 
is 

12^ .z(l,,S250m™) e(z) 

N{>S2,0,rn)^J2 -ir^V (1) 

1 = 1 ^" 



in which Vj is the comoving volume in which the j'th 
SLUGS source could have been detected in the original 
survey from which it was selected and E(z) is the ratio 
of the comoving number-density of sources of a given lu- 
minosity at a redshift z to the number at zero redshift. 
The upper limit of the integral is the maximum redshift 
at which the j'th source could be placed and just be de- 
tected above the 250-/im flux limit. Lj is the luminos- 
ity of the j'th SLUGS source, which can be estimated 



at the appropriate rest-frame wavelength at each red- 
shift (for the observed wavelength of 250 /im) from the 
two-component dust model. In the specific case of the 
submiUimet re waveband, we can u se the two-component 
dust model (jPunne fc Bales 20011 ) to estimate the lumi- 
nosity at any wavelength. If we assume that rather than 
density-evolution we have luminosity evolution, in which 
the number of galaxies stays the same but their lumi- 
nosities evolve, the equation becomes: 

Ni> ^250^™) = E / (2) 

In this case, all the information about cosmic evolution 
is included in the upper limit because 



L,{z) ^ E{z)L,{0) (3) 

in which E{z) is the ratio of the luminosity of a galaxy at 
redshift z to the luminosity of the galaxy at zero redshift. 
These equations can be modified in a straightforward 
way to estimate the source counts at any submiUimetre 
wavelength and also the cosmic background radiation. 

In practice, we have used the simple luminosity- 
evolution model from Rowan-Robinson (2001), in which 
the rest-frame monochromatic luminosity over the entire 
far-IR/submillimetre waveband is assumed to evolve in 
the following way: 

L(0=L(to)(^)"e^(^-i^ (4) 
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TABLE 2 
Predictions 



Wavelcngth/rcdshift range 


SLUGS 


Lagache 


liU /^m, total 


76944 


65973 


Tin . ^ ^ r\ ^ 
iW ^m, Z < U.l 


15489 


7267 


liU ^m, z < U.o 




iooUz 


170 /im, total 


68740 


56351 


170 fim, z < 0.1 


10088 


7049 


170 iim, z < 0.3 


26761 


17883 


250 /im, total 


222061 


170783 


250 Aim, 2 < 0.1 


12500 


13321 


250 fj.m, z < 0.3 


40073 


34887 


350 ^im, total 


50422 


32636 


350 Aim, 2 < 0.1 


2967 


3192 


350 iim, 2 < 0.3 


7389 


8170 


500 /im, total 


9734 


6998 


500 ^m, 2 < 0.1 


856 


444 


500 Aim, 2 < 0.3 


1714 


1118 



Note. — Reading from the left, the columns arc: the wave- 
length and the redshift range ('total' means all rcdshifts); the predic- 
tion of the model based on the SCUBA Local Universe and Galaxy 
Survey (see text for details); the predic tion of the Lagache model 
lILagache et al. 20031 : ILagache et al. 200l) 



in which t is the time from the big bang and to is the 

time at the current epoch. P and Q are parameters of 
the model, and we found that P = 3 and Q — 9 pro- 
duced acceptable fits to the spectral sha pe and intensity 
of the cosmic background radiation (Fi xsen et al. 199a ) 
and to the SCUBA 850 /im and Spitzer TOfim source 
counts (Coppin et al. 2006, Frayer et al. 2006a, b). 

Table 2 shows the total number of sources that this 
model and the models of Lagache et al. (2003,2004) pre- 
dict should be detected by H- ATLAS, including the num- 
ber at z < 0.1 and at z < 0.3. Figure 3 shows the redshift 
distributions predicted by the two models. The numbers 
of sources predicted by the models at the five wavelengths 
agree fairly well, although there are large differences be- 
tween the predicted redshift distributions. This differ- 
ence is not surprising because the observational data 
used to constrain the models consists almost entirely 
of number counts, which inform us about the numbers 
of sources in slices of the luminosity-redshift plane but 
nothing about the distribution of redshifts within each 
slice. Both sets of models do agree, however, in pre- 
dicting that the H- ATLAS will detect a large number of 
sources in the relatively nearby Universe (z < 0.3). 

The advantage of the models based on SLUGS is that 
it is possible to predict the properties of the H- ATLAS 
sources in any waveband in which the spectral energy 
distributions (SED) of the SLUGS galaxies have been 
measured. For example, the number of H- ATLAS galax- 
ies that are predicted to have B-band optical magnitudes 
in the range Bi <B<B2 is given by 

N{B^ < S < B2) = V / —dV 

j = l"'3™,„(Si,B2,S250f,„>,ij) 

in which the limits of the integral are the maximum 
and minimum redshifts at which the j'th galaxy would 



fall both within the optical limits and above the 250- 
/xm limit. In calculating these limits, we have assumed 
that the far-IR/submillimetre luminosity of the SLUGS 
galaxy is evolving with cosmic time in the way given by 
equation 4 but that its luminosity in the other waveband 
is not evolving. This assumption is almost certainly not 
correct, but given our lack of knowledge of the cosmic 
evolution in most wavebands it seems the safest (and sim- 
plest) one to make. It is also a conservative assumption, 
in the sense that the model probably underestimates how 
bright the H- ATLAS sources will be in the other wave- 
bands. The SED of each galaxy is necessary to relate 
the flux at the observed wavelength — in this case the B- 
band — to the luminosity of the galaxy at the wavelength 
at which the radiation was emitted. Similar equations 
can be written for any waveband. 

Figures 4-8 show the predictions for the H-ATLAS 
galaxies in five bands. The panel on the left-hand side 
shows the prediction for all the H-ATLAS galaxies and 
the one on the right for those galaxies at z < 0.3. 

Figures 4 and 5 show the predictions for the optical 
B and r bands. Several major redshift surveys have 
been based on catalogues defined in these bands and in 
§4 we use these predictions to estimate the fraction of 
the H-ATLAS galaxies that will already have spectro- 
scopic redshifts. In making the predictions for the B 
band, we have used the B-band photometry that exists 
in the NASA Extragalactic Database (NED) for most 
of the SLUGS galaxies, and we made the assumption 
that all the SLUGS galaxies have an SED typical of an 
Sbc galaxy (jColeman et al. 19801 ). Note that although 
the assumption of a single optical SED for all H-ATLAS 
galaxies is a crude one, it should introduce little error at 
low redshift, where the observed and rest frames are very 
close in wavelength. We have made the r-band predic- 
tions using the V-band photometry, which exists in NED 
for 75% of the SLUGS galaxies, with a small correction 
to the r-band using the colour transformations in Jester 
et al. (2005). 

Figure 6 shows the predictions for the K-band. We 
have used the K-band photometry in NED that exists 
for almost all the SLUGS galaxies (mostly from 2MASS) 
and again assumed a standard Sbc SED. 

Figure 7 shows the predicted 1.4-GHz continuum ra- 
dio fluxes of the H-ATLAS galaxies. We have made 
the prediction for the continuum radio fluxes using the 
1.4-GHz flux measurements that exist for almost all the 
SLUGS galaxies plus the assumption that all the SLUGS 
galaxies have a power-law radio continuum [S oc z^~") 
with a = 0.7. In this case, it seems quite likely that 
the assumption of no evolution in the radio waveband 
is incorrect because of the strong correlation between 
the far-infrared/submllimctrc and radio emission from 
star-forming galax i es both at low a nd and high redshift 
(jHelou et al. 19861 : llbar et al. 2"009l) . We have therefore 
assumed that radio luminosity also evolves in the way 
given by equation 4. 

Figure 8 shows the prediction for the HI line fluxes 
of the H-ATLAS galaxies. We used the HI line fluxes 
for the SLUGS galaxies given in NED. We have assumed 
that the amount of atomic hydrogen in a galaxy does not 
evolve, which may mean that the predicted HI line fluxes 
of the H-ATLAS galaxies are underestimates. 

The figures show that, in contrast to the sources de- 
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Fio. 4. — The histogram of B-band optical magnitude predicted for the H-ATLAS galaxies using the method described in the text. On 
the left is the prediction for all the galaxies and on the right is the prediction for those at z < 0.3. The vertical dashed line shows the 
magnitude limit of the 2dF Galaxy Redshift Survey. 



Fig. 5. — The histogram of r-band optical magnitude predicted for the H-ATLAS galaxies using the method described in the text. On 
the left is the prediction for all the galaxies and on the right is the prediction for those at z < 0.3. The vertical dashed line on the left 
shows the limiting magntiudc of the imaging part of the Sloan Digital Sky Survey; the one on the right shows the approximate limit that 
will be reached by the ESO public survey KIDS. 
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Fig. 6. — The histogram of K-band optical magnitude predicted for the H-ATLAS galaxies using the method described in the text. On 
the left is the prediction for all the galaxies and on the right is the prediction for those at z < 0.3. The left -hand vertical dashe d line 
shows the limiting magntiude of the Large Area Survey that is part of the UKIRT Infrared Deep Sky Survey 1 Warren et al. 20071) . The 
right-hand vertical line shows the limiting magnitude of the ESO public survey, VIKING. 



log,o(Flux/Jy) 



log,„{Flux/Jy) 



Fig. 7. — Histogram of the 1.4-GHz radio fluxes predicted for the H-ATLAS galaxies using the method described in the text. On the left 
is the prediction for all the galaxies and on the right is the prediction for those at z < 0.3. The vertical dashed line shows the approximate 
5(T limit of the radio survey we are carrying out with the GMRT, translated to 1.4 GHz using the assumption of a power-law radio spectrum 
(S oc u-°') with a = 0.7 (§4.3). 
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O9lo(^y km/s) 

Fig. 8. — Histogram of the 21-cm line fluxes predicted for the H-ATLAS galaxies using the method described in the text. On the left is 
the prediction for all the galaxies and on the right is the prediction for those at z < 0.3. The vertical dashed line shows the approximate 
5(7 limit of the HI survey described in §4.3. 



tcctcd in the SCUBA surveys at 850 /im, it should be 
fairly easy to follow-up the H-ATLAS sources with ob- 
servations in other wavebands. 

4. MULTI-WAVELENGTH DATA 

We selected our fields partly because of the low 'cir- 
rus emission' from dust in our own galaxy (Fig. 2) 
but mainly because the complementary multi-wavelength 
data is better than for any field of similar size. We de- 
scribe here the mult-wavelength datasets, both ones that 
exist now and ones that are likely to soon exist. 

4.1. Spectroscopy 

There have been three major recent redshift surveys: 
the Sloan Digital Sky Survey (SDSS - York et al. 2000); 
the 2dF Galaxy Redshift Survey (2dFGRS - CoUess et 
al. 2001) and the Galaxy And Mass Assembly Survey 
(GAMA - Baldry et al. 2009, Driver et al. 2009). 
The GAMA survey is still underway. All of the H- 
ATLAS fields are covered by one or more of these sur- 
veys. We have used the models from the last section 
to predict the number of H-ATLAS sources that will 
already have redshifts from one of these surveys. In 
making these predictions, we have made the assump- 
tions that (a) the 2dFGRS m easured redshifts fo r 93% of 
the galaxies with B < 19.6 (|Colless et al. 20011 ) (b) 
that the SDSS measured redshifts for 94% of galaxies 
with r < 17.77 (jStrauss et al. 20021 ) and (c) that GAMA 
will measure re dshifts for all galaxies with r < 19.4 
(jPriver et al. 20 09) . Table 3 lists the predicted numbers 

We have made the trasnformation from the bj magnitude 
system using the relationship bj = B~0.28{B — V) given in Maddox 
et al. (1990). 



of sources which will already have spectroscopic redshifts 
in the different fields. The table shows that ~ 2 x lO'' H- 
ATLAS galaxies are likely to already have spectroscopic 
redshifts, and that the percentage of H-ATLAS galaxies 
at z < 0.1 with spectroscopic redshifts is likely to be very 
high. All three surveys, of course, contain far more infor- 
mation about the galaxies than only the redshifts, such 
as measurements of line ratios, line equivalent widths, 
kinematics etc. 

4.2. Imaging from the ultraviolet to the near infrared 

At ultraviolet wavelengths, approximately half the to- 
tal survey area has been observed with GALEX. This 
has mostly been part of the GALEX Medium Imaging 
Survey (MIS), which has an exposure time of 1500s and 
a hmiting AB magnitude of ~23. We and the GAMA 
team have been awarded time for a proposal to complete 
the GALEX coverage of the GAMA fields to the MIS 
depth (P.I. Tuffs). 

In the optical waveband, the GAMA fields and the 
NGP field have been surveyed in five passbands by the 
Sloan Digital Sky Survey. The limiting magnitude of 
the r-band imaging is shown in Figure 5. As yet, the 
SGP field has only been surveyed with the UK Schmidt 
Telescope using photographic plates. In the next three 
years, the H-ATLAS fields should be covered by three 
new optical surveys. The GAMA and SGP fields will 
be observed in four passbands as part of the Kilo De- 
gree Survey (KIDS), an ESQ public survey that will be 
carried out the with VLT Survey Telescope (VST). The 
sensitivity limits of KIDS will be about 2 magnitudes 
fainter than the SDSS limits. The NGP and GAMA 
fields will be observed in five bands by Pan-STARRSl 
and the SGP and GAMA fields will be observed in six 
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TABLE 3 

Redshift Surveys 



H-ATLAS field Datascts No. at ^ < 0.1 No. at ^ < 0.3 



NGP SDSS 3366 (2895) 10865 (4563) 

GAMA GAMA,SDSS,2dFGRS 3366 (3366) 10865 (9456) 

SGP 2dFGRS 5610 (4993) 18106 (8148) 

All fields ... 12342 (11254) 39833 (22167 ) 



Note. — Reading from the left, the columns arc: the H-ATLAS fields; the rcdshifts surveys covering this field; the predicted number of sources 
at z < 0.1 in this field with the predicted number with spectroscopic redshifts in brackets; the predicted number of sources at .s < 0.3 in this field 
with the predicted number with spectroscopic redshifts in brackets 



bands by SkyMapper (jKeller et al. 20071 ). The approx- 
imate 5(7 limits for all the optical surveys are given in 
Table 4. A httle further off in time, the SGP field wiU 
be observed as part of the Dark Energy Survey. 

In the near infrared, the GAMA and NGP fields 
are being surveyed in four passbands as part of the 
Large Area Survey (LAS), a legacy survey being car- 
ried out as part of the U KIRT Infrared Deep Sky Sur- 
vey (jWarren et al. 2007[ ). The GAMA and SGP fields 
will soon be observed in five passbands as part of the 
VISTA Kilo-Degree Infrared Galaxy Survey in the In- 
frared (VIKING), an ESQ public survey that is currently 
being carried out with the Visible and Infared Telescope 
for Astronomy (VISTA). The approximate limits of both 
surveys are given in Table 4 and shown in Figure 6. 

4.3. Radio Surveys - continuum and HI 

The H-ATLAS fields have been surveyed at 1.4 GHz 
as part of the NRAO VLA Sky Survey. This survey 
is not sensitive enough, however, to detect a significant 
percentage of the H-ATLAS sources. We are therefore 
currently completing a radio survey (P.I. Jarvis) of the 
GAMA fields at 325 MHz with the Giant Metre-wave 
Radio Telescope (GMRT), which will reach an approxi- 
mate 5(7 sensitivity of 1 mJy. The limit of this survey is 
shown in Figure 7. The NGP field is also high on the list 
of targets for an early survey with the Low Frequency 
Array for Astronomy (LOFAR). 

There is currently no HI survey of any of the H-ATLAS 
fields that could yield HI measurements of a significant 
fraction of the H-ATLAS sources. However, the SGP 
and GAMA fields are natural targets for the new south- 
ern radiotelescopes that are being built as prototypes for 
the Square Kilometre Array (and for the SKA itself, of 
course). We have designed a potential HI survey of the 
SGP field using the parameters given for the Australian 
Square Kilometre Pathfinder Array (ASKAP) by John- 
ston et al. (2007). We estimate that in one month of 
observing time it should be possible to carry out a sur- 
vey of the SGP sensitive to galaxies out to z ~ 1 with an 
approximate sensitivity to HI flux of ~ 10~^ Jy km s"^ 
(Fig. 8). 

4.4. Planck and other telescopes 

The Planck Surveyor is surveying the whole sky in nine 
passbands, two of which (350 and 550 /im) are the nearly 
the same as those that will be used in H-ATLAS. The size 
(full-width-half-maximum) of the Planck beam is ~ 10 
times greater than that of H-ATLAS at the same fre- 
quency, although the sensitivity in surface brightness is 



fairly similar. Therefore, in the common area of the sky 
(one eightieth of the whole sky), the two surveys will 
be complementary, with the H-ATLAS providing high- 
resolution observations of the sources that will be de- 
tected by Planck. 

The SGP field may also be observed by the South Pole 
Telescope, a telescope designed to look for high-redshift 
clusters using the Sunyaev-Zeldovich effect. 

5. THE H-ATLAS SCIENCE PROGRAMME 

In this section we describe the six major science pro- 
grammes planned by the H-ATLAS team. We note that 
these programmes represent only a limited subset of the 
scientific projects that will be ultimately possible with 
the H-ATLAS for two reasons. First, many other projects 
will become possible as the surveys of the H-ATLAS 
fields at other wavelengths (§4) are gradually completed. 
Second, since the wavelength range from 200 to 500 /im 
is virtually unexplored and since the H-ATLAS will cover 
one eigthieth of the sky in this waveband, it is possible 
that there will be some unanticipated discoveries. 

5.1. The Local Universe 

Our models (§4) predict that the H-ATLAS wiU de- 
tect ~40,000 galaxies in the relatively nearby Universe 
(z < 0.3). Almost all the galaxies detected at z < 0.1 
and approximately half the galaxies detected at z < 0.3 
should already have spectroscopic redshifts. Most of 
these galaxies will be unresolved by the H-ATLAS beams. 
There are ~120 clusters of Abell richness class 1 or 
greater, including the Coma Cluster, in our fields. 

Apart from providing the spectroscopic redshifts nec- 
essary to calculate the intrinsic properties of many of 
the sources (luminosities, dust masses etc.), the com- 
plementary data will be helpful in other ways. Three 
are particularly important. First, the redshift surveys 
will allow us to determine the position within the cosmic 
web of each Herschel galaxy and to measure the density 
of the galaxy's environment (exactly how to do this best 
is debateable — e.g. Balogh et al. 2004). Second, the 
existence of catalogues at other wavelengths means that 
the H-ATLAS will contain more information than the in- 
dividual source detections. By coadding ('stacking') the 
Herschel emission at the positions of objects in different 
classes, we can study the dust and dust-obscured star 
formation i n objects that would be too faint to detect in- 
dividually (|Dole et al. 20061 : IDve et al. 2007bf) . An im- 
portant example is elliptical galaxies (§1). Third, we 
will be able to answer the question of what fraction of 
the optical light from a galaxy is obscured by dust. Op- 
tical astronomers have struggled with this issue for 50 
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TABLE 4 

Optical and Near Infrared Data 



H- ATLAS field Datascts uvgrizZyYJHK 



NGP 


SDSS 


22.0 




22.2 


22.2 


21.3 


20.5 














NGP 


Pan-STARRSl^ 






24.1 


23.5 


23.4 


22.4 




21.2 










NGP 


LAS 


















20.87 


20.55 


20.28 


20.13 


GAMA 


SDSS 


22.0 




22.2 


22.2 


21.3 


20.5 














GAMA 


KIDS 


24.0 




24.6 


24.4 


23.4 
















GAMA 


Pan-STARRSl^' 






24.1 


23.5 


23.4 


22.4 




21.2 










GAMA 


SkyMapper"^ 


22.9 


22.7 


22.9 


22.6 


22.0 


21.5 














GAMA 


LAS 


















20.87 


20.55 


20.28 


20.13 


GAMA 


VIKING 














23.1 




22.4 


22.2 


21.6 


21.3 


SGP 


KIDS 


24.0 




24.6 


24.4 


23.4 
















SGP 


SkyMapper'' 


22.9 


22.7 


22.9 


22.6 


22.0 


21.5 














SGP 


VIKING 














23.1 




22.4 


22.2 


21.6 


21.3 



Note. — Reading from the left, the columns are: the H-ATLAS field; the optical or near-IR survey; the quoted sensitivity limits in AB magnitudes 
(These are 5(T limits for point sources, except in the case of the SDSS in which the limits represent the 95% detection repeatability for point sources). 
^ The limits are for the planned three-year survey with Pan-STARRSl. 
^ The limits are for the planned six-epoch survey by SkyMapper. 



years (|Holmberg 19581 |D avies and Burstein 1995f ) with 
limited success, largely because optical galaxy catalogues 
have large selection effects caused by dust. The impli- 
cations for extragalactic astronomy are potentially large, 
with a recent study concluding that the optical luminos- 
ity function is significantly altered by dust extinction and 
that even bulges su ffer as much as 2 ma g of extinction at 
certain inclinations ([Driver et al. 20071 ). The ultraviolet, 
optical and near-infrared data that exist for some of the 
fields (§4) will allow us to address this issue by the sim- 
ple energy-balance technique of comparing the total dust 
emission to the total unobscurcd starlight. The only as- 
sumption in this technique is that the absorption by the 
dust is isotropic. This is almost certainly untrue but by 
averaging over many objects we can largely eliminate the 
effects of anisotropy (Driver et al. 2008). 
We give the following four projects as examples: 

• We will make the first accurate estimate of the local 
submillimetre luminosity and dust-mass functions 
down to dust masses of ~ lO**-^ M© (Fig. 9). As 
we will have Herschel measurements for ~10,000 
galaxies at z < 0.1, we will be able to compare 
the luminosity and dust-mass functions for differ- 
ent classes of galaxy (e.g. different Hubble types, 
high and low environmental density etc.). We will 
also be able to extend the study of the luminos- 
ity and dust-mass functions to higher dimensions, 
for example by examining how the space density 
of galaxies depends on both dust mass and stellar 
mass. 

• Using the results from the redshift surveys, we 
will investigate how the dust-obscured star forma- 
tion depends on the local and large-scale environ- 
ment. This comple ments prev ious optic al studies 
(|Balogh et al. 20041 : Ikauffma nn et al. 2 0041. be- 
cause our observations will be much more sensi- 
tive to starbursts, which may well have a different 
environmental dependence from quiescent star for- 
mation. 

• We will investigate how the dust content of 
the Universe and dust-obscured star formation 



has changed during the last three billion years. 
This will finally follow up an important dis- 
covery from IRAS that there is strong evolu- 
tion in the luminosity function at a suprisingly 
low redshift ([Saunders et al. 1990t ). a phenomenon 
which also may ha ve been seen in the SDSS 
([Lovedav et al. 20041 ). 

• One of the strongest correlations in astronomy 
is that between the far-infrared and nonther- 
mal radio emiss ion of galaxies (Hc lou et al. 19861 : 
llbar et al. 200 9'). A widely-accepted explanation 
is that both the far-infrared emission and the ra- 
dio emission are ultimately caused by young stars, 
with the far-infrared emission being from the dust 
heated by young stars and the radio emission be- 
ing nonthermal emission from the relativistic elec- 
trons generated in the supernova remnants pro- 
duced when the stars reach the end of their lives. 
Nevertheless, many of the properties of this rela- 
tionship are hard to explain, in particular the ap- 
proximate unity slope when the two quantities ar e 
plotted on logarithmic axes (jVlaha kis et al. 20 071). 
By investigating this relationship as a function of 
environment, redshift, star- formation activity and 
other properties, we will try to understand better 
its fundamental cause. 



5.2. Planck Point Sources and Diffuse Emission 

The Planck High- Frequency Instrument (HFI) will sur- 
vey the whole sky in six bands (3, 2.1, 1.4, 0.85, 0.55, 
0.35 mm), the first survey of the whole sky at these 
wavelengths. Apart from its measurements of the pri- 
mordial cosmic background radiation (CMB), the HFI 
will detect extended foreground emission, such as dust 
in the G alaxy, and it is also likely to detect many point 
sources ([de Zotti et al. 2 005). including nonthermal ra- 
dio sources, thermal dust emission from nearby galax- 
ies and Sunyaev-Zeldovich (SZ) sources associated with 
distant clusters, the consequence of the scattering of 
the CMB radiation by the electrons in the intracluster 
medium. 
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Fig. 9. — The triangles and stars show estimates of the dust-mass function (the space-density of galaxies as a function of dust mass) from 
Vlahakis, Dunne and Eales (2005). The vertical lines along the bottom show our estimates, on the assumption of Poisson statistics, of the 
itlcr errors on the estimate of the dust-mass function from the H-ATLAS. 



The SZ sources detected by Planck are potentially of 
great importance for cosmologists, because the number- 
density of clusters as a function of rcdsh ift depends 
critically on the cosmol ogical model (Bartel mann 20011 : 
iBartelmann and White 2002 ). It might also be possible, 
because the spectral shape of the SZ effect depends on 
the peculiar mot ion of the cluster, to measure bulk flow s 
in the universe (jKashlinskv and Atrio-Barandela 2000f ). 
which is another critical test of the cosmological 
paradigm. A major problem, however, is the contami- 
nation of the SZ effect by thermal emission from dusty 
galaxies within the large (5-10 arcmin) Planck beam, 
because there is evidence that even at moderate red- 
shift the combined emission from dust in cluster galax- 
ies is comparable to the SZ effect in the 0.85-mm band, 
and much gre ater in the two shorter wavelength bands 
(jZemcov 2007D . There is zero SZ effect at 1.4 mm, which 
leaves only the two long- wavelength bands with little con- 
tribution from dusty galaxies, but as these also have the 
worst angular resolution there is the additional problem 
of confusion with nonthermal radio sources. 

The H-ATLAS will survey one eightieth of the sky in 
the two Planck bands with the highest frequencies, but 
with much better resolution. The sensitivity to extended 
regions of low surface brightness of the two surveys is 
likely to be similar. The combination of the two surveys 
will make possible several joint projects: 

• The effects of confusion on the Planck point source 
catalogue are likely to be large, particularly in the 
highest frequency channels. These include Edding- 
ton bias ( Edding ton 1940), which can lead to both 
overestimates of the fluxes of the sources and spu- 
rious detections, and a large and uncertain effect 
on the positions of the sources. The H-ATLAS will 
provide a sample of sources with accurately known 
positions and fluxes, which can then be checked 



against the Planck point-source catalogue for the 
same area of sky for the two highest frequency 
channels. This comparison will provide an estimate 
of the correction factor for the Planck fluxes, an 
estimate of the fraction of the Planck sources that 
are spurious, and measurements of the positional 
errors of the Planck sources. 

Some of the Planck science goals, such as the detec- 
tion of t he anisotropics in the cosmic infrared back - 
ground f Lagache et al. 20031 ILagache et al. 2004[) . 
require the successful removal of galactic cirrus 
emission. The galactic cirrus should be easier to 
distinguish in the H-ATLAS images because of the 
better angular resolution, and so we will be able to 
test the success of the Planck component separa- 
tion techniques. 

The combination of the high-resolution images pro- 
vided by H-ATLAS and the spectral energy distri- 
butions of the Planck sources should make it possi- 
ble to determine unambiguously the nature of each 
source in the overlap region between the two sur- 
veys. The detailed information provided by the 
H-ATLAS for the Planck sources in this region will 
make it possible to develop better statistical meth- 
ods based on their spectral energy distributions for 
determining the nature of the Planck sources in the 
region not covered by H-ATLAS. 

In the overlap region between the two surveys, the 
H-ATLAS will make it possible to determine the 
contribution of dusty galaxies to each SZ source 
and then subtract these galaxies' emission from the 
SZ signal. The detailed information in the over- 
lap region will also provide statistical information 
about the average level of contamination for the 
Planck SZ sources not covered by H-ATLAS. 
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5.3. The H- ATLAS Lens Sample 

In principle, gravitational lensing is a powerful way of 
investigating the evolution in the mass profiles of galax- 
ies, a fundamental test of models of structure formation. 
In practice, it has proved very hard to assemble the nec- 
essary large sample of lenses. The most ambitious pro- 
gramme to date searched for lenses among flat-spectrum 
radio sources but, after high-resolution radio observa- 
tions, found only 22 le nses out of 160 00 radio sources — a 
success rate of 0.14% () Browne "2001 . 

Submillimetre surveys are possibly the ideal way 
to find lenses. At z > 1, the monochromatic flux density 
at the long-wavelength end of the submillimete wave- 
band depends on luminosity but is approximately inde- 
pendent of redshift, a consequence of the characteristic 
spectral energy distribution of dust in galaxies. A result 
of this and the strong cosmic evolution in the submilime- 
tre waveband is that sources in submillimetre surveys 
tend to lie at high redshifts (§1) and are thus likely to 
have a large optical depth to lensing. Because of the 
approximate independence of flux and redshift, a bright 
source in a submillimetre survey must necessarily have a 
high luminosity unless it is is at z < 1. Above a criti- 
cal flux density, the luminosities of any sources beyond 
this redshift must lie above the luminosity at which the 
luminosity function declines steeply — and so any sources 
above this flux density are highly likely to be lensed sys- 
tems. 

The model in Figure 10, for example, predicts that at 
Ssoo^iin > 100 mJy the sources should be a mixture of 
lensed high-redshift gal axies, nearby galaxie s and flat- 
spectrum radio sources (|Negrello et al. 2007D . Since the 
latter two catagories would be easy to remove (by using 
the submillimetre flux ratio or the presence of a bright 
galaxy), the lens yield of such a sample would be close 
to 100%. These models predict that the H-ATLAS will 
contain ^ 1500, 800 and 350 strongly-lensed galaxies at 
250, 350 and 500/im, respectively. 

This is the one of the H-ATLAS science programmes 
for which follow-up observations will be necessary, both 
to determine whether a bright 500-^m source is in- 
deed a lensed system and to scientifically exploit the H- 
ATLAS lens sample. Fortunately, a source with Ssoo/jm 
100 mJy would be be easy to map; a 30-minute observa- 
tion with the Submillimetre Array, a one-hour observa- 
tion with the VLA or a ~l-min observation with ALMA 
would be enough to confirm a source is a lens and to map 
the image structure. The prediction that large-area sub- 
millimetre surveys are the ideal way to construct large 
samples of lenses may not be correct, but if it is, there 
are many possible uses for a large sample of lenses. These 
include: 

• an investigation of the evolution of the profiles of 
the lenses (probably mostly elliptical galaxies). In 
principle, it should be poss ible to use standard 
techniques (jPve et al. 2007a[ ) to reconstruct sep- 
arately the structures of the dark-matter halo and 
baryonic component of each lens. 

• a study of the structures of high-redshift dust 
sources by reconstructing the original (unlensed) 
structure of the source (Fig. 10). Since our 
systems are lensed, we will be able to study 



galaxies well belo w the Herschel confusion limit 
()Smail et al. 19971 ). 

5.4. Active Galactic Nuclei 

The discovery that most nearby galaxies contain a 
black hole and the mass of the black hole is strongly 
correlate d with the mass of t he surrounding spheroid 
of stars (jMagorrian et al. 19981 ) was possibly one of the 
most important discoveries in extragalactic astronomy 
in recent years, because it implies the formation of the 
stars and the central black hole in a galaxy are con- 
nected. Previous submillimetre observations of high- 
redshift quasars hav e detected c±5-10% of these objects 
(jPriddev et al. 20031 iBeelen et al. 2006( 1 . and the spec- 
tral energyy distributions of the detected quasars imply 
that most of the submillimetre emission must be coming 
come from a starburst surrounding the active nucleus 
rather than directly from the nucleus itself. This sup- 
ports the idea that formation of the stars and the forma- 
tion of the central black hole in a galaxy are connected. 

We will use the H-ATLAS to detect both individual ac- 
tive galactic nuclei that fall in our fields and to carry out 
statistical analyses of AGN that are too faint to detect in- 
dividually by coadding ('stacking') the Herschel emission 
at the positions of all the AGN in a given class. Let us 
consider as an example Type 1 quasars. Figure 11 shows 
quasars drawn from the Palomar Green (PG) survey and 
the SDSS as a function of redshift and optical luminosity. 
Because of the strong correlation between luminosity and 
redshift in any flux-limited sample, it is crucial to observe 
more than one sample, so that one can compare quasars 
with the same optical luminosity at different redshifts 
and vice versa. The de tection rate of P G quasars by 
IRAS and ISO is -80% (|Haas_et al. 2003). We will use 
the H-ATLAS to extend the far-IR/submillimetre study 
of quasars to the more distant SDSS quasars, of which 
there are ~ 10'* in the H-ATLAS fields. Using the results 
of a pilot study with the Spitzer legacy survey SWIRE 
(jSerjeant et al. 2009( 1. we estimate that we will detect 
~440 quasars at z < 3 and ~ 210 quasars at z > 3. 
This is ~15 times greater than the number of existing 
detections of high-redshift quasars. We will investigate 
the properties of the quasars that are too faint to detect 
individually by a stacking analysis. 

5.5. Large-Scale Structure 

The H-ATLAS will detect -200,000 sources with a me- 
dian redshift of ~1 and will therefore contain a large 
amount of information about large-scale structure upto 
a scale of ~1000 Mpc at z '--^ 1. Many large-scale struc- 
ture projects, such as searches for baryonic oscillations, 
will only become possible in the future, as the near-IR 
(VIKING,UKIDSS) and optical (KIDS, Pan-STARRS, 
the Dark Energy Survey) surveys eventually provide pho- 
tometric redshift estimates for most of the sources. There 
are two projects that are possible immediately: 

• The H-ATLAS will make possible measurements 
of the angular correlation function on very large 
scales. These measurements will make it possible 
to discriminate between some models of galaxy for- 
mation. Figure 12, for example, shows the results 
of using a Monte-Carlo simulation to predict the 
amplitude of the angular correlation function of 
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Fig. 10 . — The left-h and panel shows the predicted number of sources as a function of 500 fim flux density plotted over the results from 
BLAST (INegrell o et al. 20071). The key to the curves is as follows: red solid — unlensed high-z protospheroids; red dashed — lensed high-z 
protospheroids; green — blazars; blue — nearby late-type galaxies. The vertical dashed line shows the approximate limit of H-ATLAS at 500 
/im. The middle panel shows a simulated image of a strongly lensed lOOmJy 500/im H-ATLAS source as observed using the Submillimetre 
Array at 870/jm in extended mode over a 6h track. The panel on the right shows the reconstructed source obtained by applying the 
reconstruction method of Dye et al. (2007a) to the simulated ring image. The image caustic of the best fit lens model is over-plotted. 
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Fig. 11.— The distribution of PG (circle) and SDSS (point) 
QSOs as a function of redshift and luminosity. 

H-ATLAS sources using three alternative models 
of galaxy formation: the hybrid model from Van 
Kampen et al. (2005) and two new models (Van 
Kampen et al. in preparation). The simulation 
shows that it should be easy to discriminate be- 
tween these models. 

• The individual sources, however, only represent 
~10% of the extragalactic background radiation 
at the Herschel wavelengths and the unresolved 
background contains a wealth of further informa- 
tion about the clustering properties of dusty galax- 
ies. A powerful technique to extract further in- 
formation is to investigate the clustering proper- 
ties of the intensity distribution on the maps after 
the high signal-to-noise sour c es have been remove d 
(jAmblard and Coorav 2007t ILagache et al. 20071 ). 
The strength of the fluctuations on large angular 
scales (linear regime) can be used to estimate the 
average mass of the dark matter halos containing 
the sources, whereas the strength of the fluctua- 
tions on small angular scales (non-linear regime) 
can be used to estimate the halo occupancy distri- 
bution (|Amblard and Coorav 2007,) . 

5.6. Galactic Dust Sources 

We will use the H-ATLAS to look for dust associ- 
ated with stars, especially debris disks and dust around 
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Fig. 12. — A histogram of the amplitude of the angular corre- 
lation function (A^fcy) of H-ATLAS sources for 100 Monte-Carlo 
realisations for each of three different models of galaxy formation: 
the hybrid model described by Van Kampen et al. (2005) and two 
new models (Van Kampen et al. in preparation). 



stars on the asymptotic giant branch. We will also 
study the distribution and spectral properties of dust 
in the interstellar medium. Although the H-ATLAS 
fields are at high galactic latitudes, in order to mini- 
mize the emission from dust within the Galaxy, IRAS 
images show that there is plenty of dust even at the 
galactic poles (Figure 13). The H-ATLAS will allow 
us to study the structure of an interstellar dust on an 
angula r scale c^lO times smaller than was possible with 
IRAS (|Mivilles-Deschenes et al. 2007D and will be pos- 
sible with Planck, while still having as good sensitivity 
to large-scale structure as the two other telescopes. The 
H-ATLAS wfll also detect the cold dust missed by IRAS. 

A more speculative idea is to look for prestellar cores 
and protostars at high latitudes. Since the dust disks in 
galaxies are thin, most of the high-latitude dust is within 
0.5 kpc. Therefore, we will be able to carry out a search 
for prestellar cores and protostars in our fields down to 
surprisingly low mass limits. On the assumption of a 
dust temperature of 20 K and a standard gas-to-dust 
ratio, we estimate that we should detect all protostars 
and prestellar cores down to a mass of ~ 0.002 Mq — well 
below the brown dwarf limit and in the Jupiter regime. 
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Fig. 13. — IRAS image at 100 fiin of a 10 x 10 dcg^ region around 
the soutli galactic pole. 

6. FURTHER TECHNICAL DETAILS AND DATA RELEASE 

The basic parameters of the survey were given in §2. In 
this section, we describe some technical issues that will 
affect the nature of the data that we eventually release 
to the community. 

An important feature of the survey is that we plan 
to observe all fields twice. The main purpose of this 
is to overcome the potential problem of low-frequency 
(1/f) noise, caused by either slow drifts in the tempera- 
tures of the detectors or other changes in the instrument 
electronics. These drifts in time can lead to spurious 
large-scale structure in the final maps, which would be a 
problem for some of our science projects. As long as the 
two sets of observations are in different directions, how- 
ever, it should be possible to use mappin g algorithms 
such as MADmap (|Cantalupo et al. 2009() to produce 
maximum-likelihood images of sky that are free of the ef- 
fects of this noise. Simulations suggest that as long as the 
scan direction of the two sets of observations is separated 
by at least 20 degrees, there shoul d be no features in th e 
final maps caused by these drifts (jWaskett et al. 20061 ). 

Achieving two sets of observations with different scan 
directions, however, is quite tricky, partly because of the 
severe constraints on the way that Herschel will observe 
the sky and partly because of the sheer size of H- ATLAS 
(11% of the time available for open-time key projects). 
Surveys will be carried out with Herschel by using the 
telescope to scan across the sky along a great circle, then 
moving the telescope a short distance in a parallel 'cross- 
leg' direction, and then moving the telescope back along 
a great circle, thus gradually building up a map of the 
sky. A complication of observing with SPIRE is that the 
SPIRE bolometers do not instantaenously fully sample 
the sky, and so to produce a fully-sampled SPIRE image 
it is necessary to carry out the scan in one of 24 possi- 
ble directions that are related to the six-fold symmetry 
of the array (there is not a similar constraint for PACS 
because it is a fully-sampled array). Other complications 
are that it is necessary to keep the Sun safely behind the 
telescope's sun-shield and that it is only possible to ro- 
tate the telescope completely around one axis. Waskett 
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Fig. 14.— The design of the surveys of the ATLAS NOP field and 
of the larger of the two SGP fields. In both cases, the purple and 
yellow lines show the areas that will be covered by individual Astro- 
nomical Observation Requests (AORs). The other colours outline 
the areas covered by complementary datasets in other wavebands. 
The red circle in the NOP figure shows the position of the Coma 
Cluster. 

et al. (2007) give a good introduction to the complica- 
tions of carrying out surveys with Herschel. The result 
of these constraints is that it is only possible to observe a 
given field during certain visibility windows, the duration 
of which depends on the position of the field. 

The easiest fields to observe are the GAM A fields, 
which have long visibility windows because they are the 
ones closest to the ecliptic plane. The Astronomical Ob- 
servation Requests (AOR), the units out of which the 
telescope's observing programme is built, of the GAMA 
fields each consist of two sets of observations with or- 
thogonal scan directions of an area of sky ~ 4 x 4 deg^ . 
These AORs, which will take approximately 16 hours to 
complete, will result in approximately square images, like 
those shown in Figure 2, although the orientation of the 
square will depend on the exact time at which the ob- 
servation is made. Thus the coverage within each square 
will be uniform but we cannot predict at the moment the 
angle that each square will make to the celestial equator, 
and thus the precise overlap with the area covered by the 
surveys in other wavebands (§4). 

The other two fields are more difficult to observe be- 
cause they have shorter visibility windows and because 
we need to map large contiguous areas of sky to make 
the best use of the complementary data in other wave- 
bands. As an example, we will consider the larger of the 
fields near the south galactic pole (SGP). To make the 
best use of the redshifts from the 2dF Galaxy Redshift 
Survey, the field has a small range of declination but a 
large range of right ascension (Fig. 2). For this field. 
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Fig. 15. — A simulation of the likely coverage of the larger of 
the fields near the south galactic pole. The range of exposure time 
shown by the colour table is a factor of four. 

the individual AORS each contain one set of observa- 
tions with a single scan direction. Figure 14 shows the 
layout of the AORs, although their exact orientation de- 
pends on exactly when the observations are done. The 
scan directions are either roughly along lines of constant 
declination, with the observations covering a rectangular 
area of ~ 3° in declination and ~ 12° in right ascension, 
or roughly along lines of constant right ascension, with 
the observations covering a rectangular area of ~ 6° de- 
grees in declination and ~ 4° in right ascension. We have 
set some constraints on when the AORs are observed to 
avoid the rectangles being rotated too far from a north- 
south line but it is not possible to make the constraints 
too severe because of the difficulty of scheduling such a 
large programme. Fig 15 shows an example of the cover- 
age we are likely to have for this field. The full range of 
the exposure time shown by the colour table is a factor 
of four, so the sensitivity in the apparent gaps is a factor 
of ~ 2 worse than in the places where we have the best 
coverage. Apart from in these gaps, every point in the 
field will have at least two sets of observations, with the 
scans in approximately orthogonal directions, and thus 
we will be able to use maximum-likelihood imaging al- 
gorithms to remove the effect of 1// noise. Figure 14 
also shows the design of the AORs for the field near the 
north galactic pole (NGP). The individual AORS again 
each contain one set of observations with a single scan 
direction. There are two sets of AORs with roughly or- 
thogonal scan directions, and each covering a rectangular 
area of ~ 9° x 4° with the longer side in the scan direc- 
tion. The orientation of the overall field, and the cover- 
age within the field, is still uncertain and will depend on 
the final telescope schedule. For all the fields, one of the 
legacy data products will be coverage maps of the fields. 

A useful byproduct of our observing strategy of two 
sets of observations of each point with an field is that we 
will be able to look for time-varying and moving objects 
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in the H-ATLAS data. Because of the scheduling diffi- 
culties, we have not placed any constraint on the time in- 
terval between two observations of the same place within 
the SGP and NGP fields, and thus it is only within the 
GAMA fields that we will have two sets of observations 
separated in time by a roughly constant interval. This 
interval is approximately 8 hours, which is well-suited for 
looking for asteroids but not, for example, for objects in 
the Edgeworth-Kuiper Belt. 

We intend to release the H-ATLAS data (maps and 
catalogues) to the community in a series of data releases. 
The first data release will be for a 4° x 4° area within 
the 9-hour GAMA field (Table 1), which wiU be observed 
during the Herschel Science Demonstration Phase (Oc- 
tober 15th-November 30th 2009). We currently expect 
to release the maps and catalogues for this region in May 
2010. The complexity of telescope scheduling means that 
we cannot yet give a firm date for the next data release. 
More information can be found on the Herschel ATLAS 
website (www.h-atlas.org). 

7. CONCLUSIONS 

The Herschel ATLAS is the largest open-time key 
project that will be carried out with the Herschel Space 
Observatory. It will survey 510 square degrees of the 
extragalactic sky, four times larger than all the other 
Herschel surveys combined, in five far-infrared and sub- 
millimetre bands. We have described the survey, the 
complementary multi-wavelength datasets that will be 
combined with the Herschel data, and the six major sci- 
ence programmes we plan to undertake. Using new mod- 
els based on a submillimetre survey of nearby galaxies, 
we have presented predictions of the properties of the 
sources that will be detected by the H-ATLAS. We intend 
to release the H-ATLAS data — maps and catalogues — to 
the astronomical community in a series of data releases. 
The first of these wiU be in May 2010. 
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